Hierarchical-Matrix Preconditioners for Parabolic Optimal Control Problems
نویسندگان
چکیده
Hierarchical (H)-matrices approximate full or sparse matrices using a hierarchical data sparse format. The corresponding H-matrix arithmetic reduces the time complexity of the approximate H-matrix operators to almost optimal while maintains certain accuracy. In this paper, we represent a scheme to solve the saddle point system arising from the control of parabolic partial differential equations by using H-matrix LUfactors as preconditioners in iterative methods. The experiment shows that the H-matrix preconditioners are effective and speed up the convergence of iterative methods.
منابع مشابه
RICE UNIVERSITY Time–Domain Decomposition Preconditioners for the Solution of Discretized Parabolic Optimal Control Problems
Time–Domain Decomposition Preconditioners for the Solution of Discretized Parabolic Optimal Control Problems
متن کاملDistributed Solution of Optimal Control Problems Governed by Parabolic Equations
We present a spatial domain decomposition (DD) method for the solution of discretized parabolic linear–quadratic optimal control problems. Our DD preconditioners are extensions of Neumann-Neumann DD methods, which have been successfully applied to the solution of single elliptic partial differential equations and of linear–quadratic optimal control problems governed by elliptic equations. We us...
متن کاملA Spatial Domain Decomposition Method for Parabolic Optimal Control Problems ?
We present a non-overlapping spatial domain decomposition method for the solution of linear–quadratic parabolic optimal control problems. The spatial domain is decomposed into non-overlapping subdomains. The original parabolic optimal control problem is decomposed into smaller problems posed on space-time cylinder subdomains with auxiliary state and adjoint variables imposed as Dirichlet bounda...
متن کاملA Comparison of Several Preconditioners in the Solution of Parabolic Systems in Three Space Dimensions Using DASPK and a High-Order Finite Element Method
We describe an algorithm for solving three-dimensional vector systems of parabolic partial diierential equations. Solutions are calculated using Galerkin's method with a high-order, piecewise polynomial hierarchical basis in space and the DAE code DASPK in time. A variant of DASPK that computes Jacobian vector products with a stored analytic Jacobian is compared with the default matrix-free app...
متن کاملMultigrid Preconditioning of Linear Systems for Interior Point Methods Applied to a Class of Box-constrained Optimal Control Problems
In this article we construct and analyze multigrid preconditioners for discretizations of operators of the form Dλ + K K, where Dλ is the multiplication with a relatively smooth function λ > 0 and K is a compact linear operator. These systems arise when applying interior point methods to the minimization problem minu 1 2 (||Ku − f || + β||u||) with box-constraints u 6 u 6 u on the controls. The...
متن کامل